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Abstract 

Starting with the valence bond solid (VBS) ground state of the ID AKLT Hamiltonian, we 



O 

make a partition of the system in 2 subsystems A and B, where A is a block of L consecutive 
spins and B is it's complement. In that setting we compute the partial transpose density matrix 
with respect to A, p TA . We obtain the spectrum of the transposed density matrix of the VBS 
pure system. Subsequently we define two disjoint blocks, A and B containing La and Lb spins 
respectively, separated by L sites. Tracing away the spins which do not belong to AUB, we find an 
expression for the reduced density matrix of the A and B blocks p(A, B). With this expression (in 

O" 1 

the thermodinamic limit), we compute the entanglement spectrum and other several entanglement 
measures, as the purity P = tr(p(A, B) 2 ), the negativity J\f, and the mutual entropy. 
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Entanglement is a fundamental measure of how much quantum effects we can observe 
and use, and it is the primary resource in quantum computation and quantum information 
processing [TJ [2]. Also entanglement plays a role in quantum phase transitions [2J H], and 
even it has been experimentally demonstrated that entanglement may affect macroscopic 
properties of solids [HI E]. Currently there is considerable interest in quantifying entangle- 
ment in various quantum systems. Measures of entanglement, like entanglement entropy of 
correlated electrons, fermions in conformal field theory, spin chains, interacting bosons and 
other models have been studied [THT2"] . For a review of entanglement entropy as an area law 
see [13J. 

An important measure of entanglement is negativity, introduced in [Tj5]. Negativity is 
a useful quantity to characterize quantum effects in mixed systems, where the standard 
mutual information entropy fails to provide a clear separation between classical and quantum 
correlations. Negativity is also useful in the context of quantum information because it does 
not change under local manipulations of the system [24] . It is computed from the partial 
transpose density matrix with respect to a subsystem A, p TA and essentially measures the 
degree of which p TA fails to be positive. Despite the usefulness of negativity, this quantity 
is usually difficult to compute analytically. 

In this paper, we calculate the negativity of blocks in the ground state of the spin chain 
introduced by Affleck, Kennedy, Lieb, and Tasaki (AKLT model) [14J. This state is known as 
Valence-Bond Solid (VBS). The AKLT model plays an important role in condensed matter 
physics, being the first rigorous example of an isotropic spin chain which agrees with the 
Haldane conjecture [15], i.e. Haldane's suggestion that an anti-ferromagnetic Hamiltonian 
describing half-integer spins is gapless, while for integer spins it has a gap [IB] . AKLT is also 
central in an specific scheme of quantum computation, namely measurement based quantum 
computation [T71 125] . 

While this ground state is fourfold degenerate for an open boundary chain of spin 1 at 
each site, it becomes unique for a chain consisting of bulk spin-Is and two spin-l/2's at the 
boundary [20] . 

An implementation of AKLT in optical lattices was proposed in [18] , and the use of AKLT 
model for universal quantum computation was discussed in [17] and in [26]. VBS is also 
closely related to Laughlin ansatz |2lj and to fractional quantum Hall effect |22j . 

In the first section we quickly review the formulation of the AKLT model and the VBS 
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ground state with its extension to make it unique. We also introduce the definition of the 
density matrix associated with the VBS ground state. In the second section we discuss the 
case of a bipartition of the pure ground state. We re-derive the spectrum of the partial 
density matrix pa = trgp using our simpler approach obtaining the results already shown 
in [29]. We also computed the transposed density matrix p TA to illustrate our method. For 
this case we compute the full spectrum, along with the eigenvectors of p TA . We also give 
a value for the negativity in this case, which decays to a constant value twice as fast as 
expected from the correlation functions. In the third section we define two blocks A and B, 
separated by L sites. We compute the density matrix of the mixed system AU B p(A, B), 
evaluated by tracing out the environmental degrees of freedom. We obtain the spectrum of 
p(A, B) and the entanglement spectrum as function of the separation L between the blocks 
and the size of A and B. The purity of this system corresponds to the one encountered for 
maximally mixed states (up to second order corrections). In this section we find that the 
negativity for this system vanish for non adjacent blocks. We also study the case of periodic 
boundary conditions. In the fourth section we obtain the mutual entropy of the system, in 
the limit of infinity blocks A and B. 

I. THE AKLT MODEL AND THE VBS STATE 

The one dimensional AKLT model that we will consider consists of a chain of TV spin- Is 
in the bulk, and two spin-1/2 on the boundary. The location where the spins sit are called 
sites. We shall denote by Sk the vector of spin-1 operators and by Sf, spin— 1/2 operators, 
where b — 0, N+ 1. The Hamiltonian is H = i^Buik + no,i + IL^at+I; where the Hamiltonian 
corresponding to the bulk is given by 



and the sum runs over the lattice sites. The boundary terms II describe interaction of a spin 
1/2 and spin 1. Each term is a projector on a state with spin 3/2: 




N-l 



(1) 



i=l 



n ,i = -(i + s ■ St), u N>N+1 = -(i + s N ■ s N+ i). 



(2) 
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In order to construct the ground state |VBS) of H we can associate two spin 1/2 variables 
at each lattice site and create the spin 1 state symmetrizing them. To prevent the formation 
of spin 2, we antisymmetrize states between different neighbor lattice sites. Doing this we 
are sure that this configuration is actually an eigenstate of the Hamiltonian, with eigenvalue 
(i.e. the projection of \VBS) on the subspace of spin 2-states is zero). Noting that the 
Hamiltonian H is positive definite, then we know that this is the ground state. 

We can associate a graph to this state, defining dots as spins 1/2, links as anti- 
symmetrization, and circles as symmetrization. The graph representation of the VBS 
ground state is then given by 

1 2 N N+1 

FIG. 1. Graphic representation of the ID VBS state. 



It is possible to write down a compact expression for this VBS state using bosonic 
variables. Following [22] . we make use of the Schwinger boson representation for SU(2) 
algebra at each site j, namely = djbj, Sj^ = ajb^, Sj = \{a\a-j — bjbj), with 
[S?,Sp] = ±Sfdij, [S^,S~] = +2Sf8ij, where a and b are two sets of bosonic creation 
operators, with the usual commutation relations [aj,aj] = [bi, &t] = 5ij, [ai,aj] = [bi,bj] = 
and correspondingly for and tf. This two sets commute in each and every lattice site, i.e. 
[aj,6j] = [a\, bj] = 0. To have a finite dimensional representation of SU(2), it is necessary 
to impose one more condition on a and b, given by ^(ajaj + b^bj) = Sj, with Sj the value of 
the spin at site j (in this paper we have Sj = 1 for j = 1..N and 1/2 at j = 0, N + 1). In 
this language the VBS ground state is given by [22] 

N 

|VBS}=nK t Ci-al + i&l)|0). (3) 

where |0) = <S> sites \Qa, j) <S> \0b,j)- The state \0 a ,j) is defined by a,j\0 a ,j) = 0, and it's 
called the vacuum state for the set of operators a. \0b,j) is defined similarly for the set b. 
In [20] the authors prove that this ground state is unique for the Hamiltonian H, then we 
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can construct the density matrix of the (pure) ground state p = ^^s^vbs) • This is a one 
dimensional projector on the |VBS) ground state of the Hamiltonian. 



II. DENSITY MATRIX PURE STATE 



In order to compute the partial transposed density matrix of the VBS system, we define 
three subsystems A, B\ and B2, where A is a block of L spins 1 and B = B\ U B2 is it's 
complement (see Fig. [2]) 



1 



<6K3K> 



m 



m+L-1 



<GK> 



B: 



A=L sites 



N N+1 



B: 



FIG. 2. Partition of the ID chain in three subsystems A,Bi,F>2- Subsystem A, consisting of L 
sites. Subsystem B = B\ U B2 is the complement of A. 



The partition is defined by (with 1 < L,m < N) B\ = {sites i, < % < m — 1}, 
A = {sites i, m < i < m + L — 1}, B 2 = {sites i, m + L<i<N + l} 

We can split the expression ^ in the corresponding states of the subsystems 

I VBS) = (oLA - - ct K+x tf K )\A,B). (4) 

where K = m + L — 1. \A) and \B) are the VBS states of the A and B subsystems, defined 
by \A) = \ m K ) , \B) = |o m_1 )U'+i iV+1 ) ; where the states of the form |/ J ) are defined as 

j -i 

\i J ) = X\{^\ + i-4. + Am, (5) 
1=1 

respectively. This states describe spins 1 in the bulk (i.e at I ^ I, J) and spin 1/2 in 
the boundary. To recover the spin 1 at those boundaries sites, we introduce the following 
notation (i^lY = a[, and (iplY — b\., to have 
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(rJ(^\A) = \ c A d ), {r m -i)\^K + i) ] \B) = \ c B d ) M = l,2), 

then eq. Q becomes |VBS) = l 2 ^ 1 , 1 B 2 ) - \ l A 12 B 2 ) - \ 2 A 2 , 1 B 1 ) + ^A 22 B l ). 

The four states \ a A v ) belong to the Hilbert space of the block A of length L. They span 
the kernel of H-Buik(A) [14J, but they are not orthogonal to each other. We make use of the 



classical variable method, introduced in [22] (see appendix VIA), to prove 



llVril 2 = \y- i~ l ) a+V (4) J ' < M T^> = -\ [-I)" for v?u, (6) 

valid for L > 1 {L G N). The norm ||m|| 2 is defined as usual ||"u|| 2 = (u\u). All other 
combinations vanish. We can perform a rotation of this basis in order to make the overlap 
(JgJ) vanish. The new basis is defined by 



\A ) = ^(I'A 1 ) + \ 2 A 2 )), \A t ) = ^A 2 ) + \*A% 
\M) = ^d'A 2 } - IM 1 )). \A ) = -^(l 1 ^) + \ 2 A 2 }), 



In this basis (here //, v = 0..3) the norm is given by 



11 ( P ^ 



where s M = (—1, —1,3, —1). This four different eigenstates of the bulk Hamiltonian corre- 
sponding to the block A, can be labeled by the Bell pair that is formed between the spins 1/2 
at the boundary. The boundary operators a, b which act on the subspace B, also organize 
themselves in irreducible representations, with the only condition that adjacent boundary 
operators acting on A and B cannot create a state of spin 2, as required by the VBS ground 
state symmetry. We define the following operators for further simplicity (here, sum over 
dummy variables c = 1, 2 and d = 1, 2 is assumed) 

T^^j)=^Mc d ^\ (A* = 0..3), (8) 
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with (Tp = (il, a±, 02, 03) being I the 2x2 identity matrix and Oi are the Pauli matrices. 

The operators keep explicit the symmetry between the operators a and b, remaining 
unchanged (up to phase factors) when we interchange the operators a) and . This operation 
corresponds to take the transpose of o^, so (cr^) T = o^, for fj, — 0, 1, 3 and (cx 2 ) T = — cr 2 . As 
the set of operators T M acting on the outer edges of an state is just a linear combination of the 
states defined in ([7]), the states of the form T^(i,i) \ j) are a basis for the four dimensional 
space of degenerate ground states of the bulk Hamiltonian Q. 



With the introduction of this operators, we can write the identity (see VI B ) (sum con- 
vention, with p = 0..3) 

T 2 f (z, i + l)T^(j,j + 1) = ~\V{i + l,j)Tj(i,j + 1), (9) 

Using this identity in eq. Q, we find that the VBS state has the decomposition 
(T u \m,K)\A) = \A u )) 

|VBS) = -ir M t( m - 1,^ + 1)1^,5). 
Now we can write the density matrix for the VBS state p = ^^^s|vbs^ 

p = T^(m-l,K+ l)\A^ B)(A a ,B\T a (m - 1, K + 1). (10) 
We can define the state \s) = T^(m — 1, K + 1)|^4 M , B). In terms of this state, the density 



matrix (10) takes the form p = \s)(s\, and the eigenvector is clearly \s) with eigenvalue 1. 
this is natural because so far we have just taken another basis to represent p, which was 
already a projector onto the VBS ground state. 



It's important to note that if we trace the B block in expression (10), we get the partial 
density matrix respect to the A subsystem pa = |A At )(A M |. This expression was already 
found in [29]. The von Neumman or entanglement entropy is given by Sa = — Tr(p^ In pa) = 
—\q(L)\\i\q(L) + 3Ai(L) In Ai(L) and scales with the length of the boundary as expected. 
The entanglement spectrum for pa = \A^)(A^\ is £1 = In { tj^^=l \ no degeneracy and 
£2 = In (y i_^)-L j with triple degeneracy (L 7^ 0) [29] . 

We transpose the elements belonging to the A subspace to obtain the partial transposed 
matrix p T - A 



p TA =T l J(m-l,K+ 1)| A,, B){A^ B\T a (m -1,K + 1) (11) 

where it is understood that all the operators are evaluated at the boundary sites m — 1 
and if + 1 and L — K + 1 — m > 1. With this explicit form of p TA , we can compute the 
eigenvectors and eigenvalues. 

Eigenvectors and negativity of p TA 

Defining 

VL)^(l + 3(-i) L ) and A,(L) , 1 (l - 

we have: 

Eigenvectors of p TA corresponding to Ai(L) (6-fold degeneracy) (here all the operators 
act in the sites of B that are nearest neighbors of the block A, namely = Tj(m— 1, K + 1)) 

|e 1 )=T t|A , J B); |e 4 ) = T^IA, -B) - T f \A 3 , B), 
\e 2 ) = T^\A 3 , 5); |e 5 ) = T^A, B) - 7^1 A, 5), 
|e 3 ) = T^IA, £?>; |e 6 ) = 7^1 A3, 5) + T^A, B). (12) 

Eigenvectors of p Tjv corresponding to — Ai(L) (3-fold degeneracy) 

\e 7 )=T^\A 3 ,B)+T 3 ^\A ,B), 
\e 8 )=T ^\A 1 ,B)+TS\A ,B), 

\e 9 )=Tj\A 3 ,B)-Tj\A 1 ,B). (13) 
Eigenvectors of p TA corresponding to a/A (£)Ai(L) (3-fold degeneracy) 



\e w )=Tj\A 2 ,B)-Tj\A 1 ,B), 
\e 11 )=Tj\A Q ,B)+Tj\A 2 ,B), 

\e 12 )=T^\A 2 ,B)-Tj\A 3 ,B). (14) 
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Eigenvectors of p TA corresponding to — ^X (L)Xi(L) (3-fold degeneracy) 



\e 13 )=T 1 *\A 2 ,B)+Tj\A 1 ,B), 
\e u ) = T 2 *\A <i ,B)-Tj\A 2 ,B), 

\e 15 )=Tj\A 3 ,B)+Tj\A 2 ,B). (15) 

Finally the eigenvector of p TA corresponding to Xo(L) is |eie) = T 2 ^|yl2, B). The negativity 
of the partial transposed matrix p TA obtained from the pure VBS state as the sum of negative 
eigenvalues, namely 

* = I (l - ("3r L + V(l-(-3)- L )(l + 3(-3)-^)) , 

= 3(Ai(L) + V%(L)A (£)), ( 16 ) 

valid for L > 1. This expression simplifies in the case L = 1 toAi=i = |(1 — (— 3) _1 ) = 1. In 
the thermodynamic limit [L — > oo) the negativity approach a constant value exponentially 
fast 

■^ao = ^( l-^(-fT)- (17) 

Special case L = 0. 



If we evaluate the expression (16) in the case L = 0, we get zero. This is obvious because 
in this treatment, L = means that there is no subspace A to transpose. Nevertheless there 
is another case which we have not considered so far and is analogous to take the block A 
to be empty. Namely instead of making a partition in three subspaces, we can consider a 
partition into 2 subspaces. Let's call this subspaces B% and B 2 . The site where we make the 
partition is m (with m inside the chain). Then the partition is given by 



B\ = {sites i, < i < m} B 2 = {sites i,m + l<i<N + l} 
The density matrix is trivially written as (using the notation introduced in (|6 



p=(5 a c 5 b -5 a d 5 b c )\B a 1 b B 2 ){B c 1: d B 2 



(18) 



The transposed density matrix with respect to B 2 is 

p T B2 = (5«8 d - 8 a d 8 b c )\B a 1 , d B 2 )(B c 1 b B 2 \. (19) 

The eigenvectors of p Ts 2 are, for Aj ^ = | (three fold degeneracy), \ek) = Tl(m } m + 
1)\B 1 ,B 2 ) (k = 1,2,3), while for the negative eigenvalue X 2 = — | (no degeneracy) |e ) = 
T(|(m, m + 1) l-Bi, B 2 ). The negativity in this case is Af = 1/2. 

III. DENSITY MATRIX FOR THE MIXED SYSTEM OF 2 DISJOINT BLOCKS 
A. Open boundary conditions 

So far we have studied the density matrix for the pure VBS system. In this section we 
want to extend our results to the case of mixed systems. We will study the mixed system 
composed of two blocks A and B of length La and Lb, obtained by tracing away the lattice 
sites which do not belong to these blocks in the VBS ground state. This situation is described 
in Fig [3j 

1 ... L -l L ... m-1 m ... K K+l ... M M + l ... N N + l 

v v ' v v * v v * v v * v v ' 

C (l_c sites) A (L A sites) D(L sites) B (L B sites) E (L E sites) 

FIG. 3. We made a partition of the VBS state in 5 sectors, labeled A, B, C, D and E as shown in 
the figure. To obtain the density matrix for the blocks A and B, we trace away the spin variables 
at the sites inside C, D and E. 



To define the blocks, we partition the V + 2 sites of the chain into five different subsets, 
A, B, C, D and E, of different length. 
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Given I, J, K, M, N five positive integers ordered as < I < J < K < M < N + 1, we 
define: 



• 


Block C = 


{sites 


i, < % < I — 1}, with length Lc = 


I. 


• 


Block A = 


{sites 


i, I < i < J — 1}? with length La = 


-J-I, 


• 


Block D = 


{sites 


i, J < i < K — 1}, with length L = 


K-J, 


• 


Block B = 


{sites 


i, K < i < M - 1} with length L B 


= M - K and 


• 


Block E = 


{sites 


i, M < i < N + 1} with length L E 


= N + 2-M. 



We are interested in the density matrix for the mixed system of two different blocks. In or- 
der to compute the density matrix we have to trace away the sites outside the corresponding 
blocks. To do that we use the following results for the bulk states 



(D^\D U ) = S^X^L) with \ ll =l(l + z(L)s lt ), (20) 
z{L)= {-\f- 8„ = (-1,-1,3,-1), 
(C, E]T^(I - 1, M)T V \I ~ 1, M) [C, E) = 5^. (21) 

Density matrix of the blocks C & E 

The simplest case occur when we trace away the A, D and B blocks. The density matrix 
for the C and E blocks is 



( |VBS)(VBSn (22] 
Pce - Tr ABD | (VBS|VBS) | , (22) 

= X,Tl(I - 1, M)\C, E) (C, E\T,(I - 1, M), 

= X^\[C,El)([C,El\. 

the partial transposed density matrix respect to the E system is 



P T C E E 



(A«-(-l)' 



, A2 — Ai 



)\[C,El)([C,E], 



(23) 



This is a sum of projector operators (a consequence of (20)). The negativity is non 
vanishing just for the case La + L B + L = 0, when we have Af = 1/2. 
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Density Matrix of the blocks A &; B, case L > 1. 

In this case we compute the density matrix for blocks A and B. We obtain this density 
matrix by tracing away the states on the C, D and E subspaces. 

( |VBS)(VBSn (24) 

~ ^ CDE { (VBS|VBS) J ■ (24) 



Using the identity ( |64| ) (See VI B), we can write the VBS state as a linear combination of 
products between the different four fold degenerate ground states of the bulk Hamiltonian 
([TJ) in the form: 



|VBS) = M ta/pa T^(I- l,M)\C : A p ,D u , B P ,E), with 

AW = (-1) W + - w <f w + <T W). (25) 

here we have introduced three types of tensors, the Kronecker delta symbol in 4 dimen- 
sions S a 0, the diagonal tensor g^ v = (— 1, +1, +1, +1) and the Levi Civita symbol in four 
dimensions e pupa , which is a totally antisymmetric tensor, with e pvpa = sign of permutation 
(/x, is, p, a) if (/!, v, p, a) is a permutation of (0, 1, 2, 3), and zero otherwise. 

Using this representation of the VBS state, it's easy to write down the density matrix 



(24) using the orthogonality of the bulk ground states (20). We find that the density matrix 



Pab is 



Pab = M pupa M auP(J \A p , B p ,)(A a , B p \, (26) 

with the tensor M pupa M au p a given explicitly by (summation over dummy variables v and a 
is assumed) 

/ S p — S \ 

P \ m , (27) 



with S pa = (s p + s a )/2. We can identify two parts in (26), the first term which does not 
depend on z and the rest which is linear in z. The first term is a projector on the ground 
states of the bulk of A and B, namely 

Po (A,B)=5 p J p p\A p ,B p ,)(A a ,B /3 \, (28) 
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while all the other terms, proportional to z(L), have vanishing trace. From this expression, 
we can compute the entanglement spectrum associated with p(A, B), in the limit L, Li,L 2 ^> 
1. This density matrix has rank 16, and is exponentially close to a maximally mixed state, 
the eigenvalues are (using x\ = (— 3) _Ll , x 2 = (— 3)~ L2 and z = (— 3)~ L ) 



, n 1-xx — x 2 — z 14 l + 3xi + 3x 2 + 32 

1 + Xi + x 2 + z 1 



Ai5,i6 = * 2 ± -a/^ 2 + (^1 + a>0Oi + x 2 - z). (29) 

lb o 

The entanglement spectrum of p(A, B) is ^ = — In Aj. Explicitly we have 



= 4 In 2 + X! + x 2 + z, {£ }£i 2 = 4 In 2 - 3x x - 3x 2 - 3z, 
65,16 = 4 In 2 - xi - x 2 - z + \J z 2 + (^l + ^2X^1 + x 2 - i). 



(30) 



The purity, defined as 7 = Tr(p 2 ) corresponds in this limit to the purity of a maximally 
mixed state, up to terms of second order in xi,x 2 and z. We have 7 = i + 0(2). The 
general results for arbitrary L, Li, L 2 are given in appendix |VI C 



If we call px(A, B) to all the linear terms in z(L) on (26), we can write for brevity 



Pab = Po{A, B) + z{L) Pl {A, B). 



(31) 



From the expressions (26) and (27) we can obtain the partial transposed density matrix 



with respect to the A subsystem. 



n TA 
Pab 



bfiabpP — z{L) I — dpuSaplSpa — £p,/3ap 



\A^,B p ,)(A a ,B02) 



from where, comparing with equations (26) and (27), we learn that 



Pab( z ) = Pab(-z). (33) 

Given this result, and the fact that Pab{~ z) is also a density matrix (proved in the 
following theorem), the negativity vanishes for L > 0. 

Theorem 1. The negativity of the transposed density matrix p^(z(L)) is strictly zero for 
two blocks separated by L > 0. 
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Proof. Consider the family of density matrices pab(z) = Po(A,B) + z(L)pi(A, B), denned 



in eq. (31). Recalling that the space of density matrices is convex [27], meaning that for 
two density matrices pi, p%, the operator p = Xp 1 + (1 — X)p 2 is also a density matrix for 
A G [0,1], we proceed as follows. We take the first two members of the family pab(z), 
namely pab(zi) and pab(z 2 ) for fixed z\ = l,z 2 = —1/3 (We can take any pair different z\ 
and z 2 , but the greater z is achieved for z\ = 1, z 2 = — 1/3 (or vice versa)). By the convexity 
of the space of density matrices, p = Xpab( z i) + (1 — ^)pab{z 2 ) is also a density matrix. 



Using (31), we write explicitly p = p (A,B) + (\z\ + (1 — A)z 2 )pi(A, B). We can choose 



A = |(1 - 3(-|) L ) G [0, 1] for L > 1. Using this A, we find 

p = p (A,B)-z(L)p l (A,B). (34) 



Then p = Pab{~ z ) is also a density matrix, for L > 1. Now, by (33), p T A A B (z) is also 
density matrix for z < 1 (L > 0). Then the negativity (sum of negative eigenvalues) of 
Pab( z ) vanishes. □ 

Special case L = 0. 

As in the previous section, we analyze separately the case L = 0. In this case the block 
D is not present and we cannot take a trace over it. 



We can study this scenario using the following identity (see VI B ) 



T 2 f (2, i + l)T 2 (j,j + l)Tl(k, k + l) = (-lym^xTjii + hj)T v \j + 1, k)T x \i, k + 1) 

where m^ A = =^(g^5 2X + g^S^x ~ dxvh^ + 5^%q2a) (35) 

The VBS state splits into 



|VBS) = 7%(i, i + 1)TJ(j, j + l)T 2 T (fc, k + 1)|C, A B, E) 



(36) 



with i — Lq — 1, j = Li + Lq — 1 and k = Lq + L\ + L 2 — 1 and Lq being the number 
of sites in C, L\ the number of sites in A, L 2 the number of sites in B. Here the states 
|C), \A), \B), \E) are defined as in the previous section taking m — 1 = K (L = 0). 



Using the identity (35), the equation (36) becomes 
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|VBS) = (g^S 2 x + 9 2u Spx - gxvh» + g ua e^ 2 x)Tx\i, k + 1)|C, A M , B„, £>) 



Now, given that (see appendix VI A ) 



(C, £>|T A (i, k + l)T x S(i, fc + 1)|C, £> = 5 AA ,, (37) 
we can write the normalized density matrix Tr c,eP = Pas as 

PAB = \&p,a&pp ~ [finp$aP ~ fipa8[M/3}Spa + £«^p ( ~ ^ J I^M> B v )(A a , Bp\ 



This expression is analogous to the (26), with z(L) — 1 (L — 0). From this result, the 



transposed density matrix respect to A is now given by 



PP n " )\A^B v ){A a ,Bp\ (38) 

Using this expression it is possible to compute the eigenvalues of the transposed density 
matrix and the negativity. For the results at finite size of the block A (L x ) and block B 



(L 2 ) see appendix (section VIC). The negativity in the asymptotic limit L\ — > oo,L 2 — > oo 
is given by 

= )■ (30) 



2 4 \ V 3/ V 3 



B. Periodic Boundary Conditions 

Using the same technology developed previously, we can also analyze the case of periodic 
boundary conditions. This state is unique, given that the coordination number for each spin 
is two [20]. 

In this state, we make a partition in four sectors, labeled by their length as La, Lb, Lc,Ld, 
with La + L B + L c + L D = L the total length of the system. We trace away states from 
sectors that do not belong to A U B (See fig. 111). 
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FIG. 4. We trace blocks C and D, leaving a reduced density matrix in terms of the states of the A 
and B blocks. 



We split the |VBS) as in the previous section, but now the difference is that the norm of 
both states that we trace out, namely the block C and D, is nontrivial, each one contributing 
with a factor (C a \C'g) = \a(Lc)S a /3 and (D a \Dp) = \ a {.Lr))8 a p. 

The VBS state can be rewritten as |VBS) = M pupa \C a , A p , D„, B p ), with M pvpa = 
(—l) u (g^ u 5 pa + g vp 5 pa — g ua S pp + g va t papa ). The reduced density matrix in this case is 

P ab = M pupa MaupaK{.L D )\ u (L c )\A ll ,B p ){A a ,Bp\. (40) 
The tensor W ppa p = M tlvp(T M av p a \ ]/ {L D )X a (Lc) is given explicitly by 

W ppa j3 = 5 pa S p /3A a i3(Zd, Z c ) — 5 ap 5 p pT ap (Zd, Z c ) + e a /3ij, p Rpfi/3a(Zd, Z c ) — S^pSapT^a^ — Zd, —z c ), 

where z c = z{Lc) and Zd = z(Ld). The tensors A a p(x, y), T aa i(x, y) and R a 0a'i3'(%,y) are 
respectively given by 



1 + (s a s/3 + s a + se)xy s a + s a > ( x + y 

and Rawix, y) = Sa ^ ip-y). (41) 



As with the case studied in section III, the partial transposed operator p T A A B is exactly 
Pab(— Zd, —z c )- With this result, a vanishing negativity of the system is analogous to prove 
that PAB^—Zd, —z c ) is a density matrix. We have 
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Theorem 2. The negativity of the transposed density matrix for the system with periodic 
boundary conditions p^(zd, z c ) is strictly zero for L c and L D ^ 0. 

Proof. Again we proceed as before. The density matrix p^# defines a family of operators 
pAB(z d , z c ) = Po + z dPl + z c p 2 + z d z c p 3 , with 



Po 

Pi 

P2 
P3 



\A„B p )(A a ,B 



1 + 3z(L) 

Sa _ S ^ f A A A A \ -i- ( g a — Sp + Sfi — s p) 
2 + 6z(L) ^ p ^ ~ °w°<x0) H TTToTTm e <*Pw 

Sa Sp, 



2 + 6z(L) 

SgS/3 + S a + S/3 

l + 3z(L) 



4 + 12z(L) 

(X X X X \ ( Sa Sfj -\- Sp S p ) 

(<W<W - W 4 +m(L) 6 

{bapbpP + fippSaP 



a/3p,p 



B p )(A a , Bp\, 
\An, B p ){A a , Bp\, 



[Sg + Sp) 

1 + 3z(L) 



At> B p )(A a: B^ 



We choose four different members of this family, p a = pab(1> 1), Pb = Pab{—\, 1), Pc = 
Pa_b(1, -3) an d Pd = Pab(— h, —3)- Recalling that the space of density matrices is convex 



we have that p = ap a + /3p b + 7p c + (1 — a — f3 — "f)pd is also a density matrix for 
< a, fi, 7 < 1. Choosing a = ^ + ^O^d - £ c - z d ), /3 = Jj + ^(9z d - 15z c + 9z c z d ) and 



7 



+ ^(9z c — 15zd — 9z c z d ), we have 



32 1 32 



p = p - z d pi - z c p 2 + z d z c p 3 = Pab(-z c , -Zd), 



(42) 



is also a density matrix for < a, /3, 7 < 1. This condition breaks down when z c or are 
equal to 1. 

□ 



IV. MUTUAL ENTROPY 

Having the explicit expression for the density matrix of two blocks and the proof of 
vanishing negativity for any separation of the blocks greater than zero, a natural question 
to ask is, Does the mutual entropy vanish in this case, too?. We found surprisingly that 
the answer is negative, i.e. there is non zero mutual entropy even when the separation L is 
greater than zero. The mutual entropy decays exponentially, as expected from the spin-spin 
correlations. 
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From (26), we can in principle compute the eigenvalues and eigenvectors of pab, where 
Pab can be written as a 16 x 16 matrix. This dimension is fixed by the dimension of the 
ground state space for each block. Being V the Hilbert space spanned by the vector \A^, B v ) 
we have 



Dim(V) = Dim(KerA) x Dim(KeiB) = 16 (43) 

where KerA is the Kernel of the bulk Hamiltonian defined on the Hilbert space of the 
block A, and similarly for Keri?. As we have shown, this spaces are spanned by the states 
{Afj), \B V ) //, v = 0..3, states which are orthogonal but not normalized in our convention. In 
the thermodynamical limit, when the length of each block goes to infinity L\ — > oo,L 2 — > oo, 
the states \A^) become orthonormal 

(A U \A,) = ^, (B,\B V ) = 5^. (44) 
In this limit, the matrix elements of pab are 



B„\p A B\A a , Bp) = S^aSvp + z^lS^Sap - 5 va 5 '^a + z(L) 
The eigenvalues of this matrix are 

Xi(z) = — (1 + 3z), 4-fold degeneracy (45) 

\n(z) = — (1 — z), 12-fold degeneracy. (46) 
16 

Using the spectral theorem, we find that the entropy of the system described by Pab, i-e., 
the entropy of two blocks of infinite length separated by L sites is 

S[A, B] = -Tt( P ab Hpab)) = 2 In 2 - |(1 - z) In (i=*) - 1(1 + 3z) In (^) . (47) 
The mutual entropy/information is defined as usual 



I(A,B) = S[A) + S[B)-S[A,B], 
18 



(48) 



the entropy of a block of length L in the AKLT model was calculated in [29] , and also can 
be obtained trivially from our results of section [TTJ In the limit of infinite length, we have 



S[A] = S[B] = 2 In 2. (49) 
The mutual information is finally 

I (A, B) = -(l-z) ln(l - z) + -(1 + 3z) ln(l + 3z), (50) 
where z was defined before as z = z(L) = (— 

V. DISCUSSION AND CONCLUSION 

In this paper we have derived the entanglement spectrum of the density matrix of two 
blocks belonging to the ID VBS state, corresponding to the ground state of the spin 1 AKLT 
Hamiltonian. The eigenvalues of the density matrix decay exponentially with the length of 
the blocks and their separation to the eigenvalues of a completely mixed state. This decay 
was expected from the behavior of the correlation functions (SqS j l ) = | (— Sy. The novel 
result is that in the thermodynamic limit, the density matrix p(A, B) eq. ([26)) is maximally 



mixed, with all the eigenvalues A, = ^ for i — 1.-16. The density matrix (26) in this limit is 
a projector on the different ground states obtained as a tensor products of the four ground 
states of A and B blocks \A,j)\B v ). 



The density matrix for this system (26) is clearly non separable. This can be rigorously 
proved, using that an state is separable iff the quantity Tr (Opab) — for any Hermitian 
operator O satisfying Tr(OP <S> Q) > 0, where P and Q are projections acting on the 
Hilbert spaces associated to subsystems A {T-La) and B (Hs)- If we choose for example 
O = r(\A ,B Q )(A 1 ,Bi \ + \A U B 1 )(A , B \), it is easy to see that Tr(0P ® Q) = 0, while 
Tr(Op) ex r, then choosing r properly we can make Tr(Op) < 0, proving that the state is 
not separable. 

The fact that the operator p(A, B) describe a state which is non separable, together with 
the result that the negativity vanishes for any separation of the blocks A and B tell us that 
this state is a bound entangled state [31] . In [31] the authors show that this kind of states can 
not be distilled by local action to create an useful entanglement for quantum communication 
tasks such as teleportation. However, bound entanglement is still of interest as it can be 
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used to generate a secret quantum key [30], as well as be used to enhance the fidelity of 
conclusive teleportation using another state [32] ■ In the context of many-body systems, it 
has so far been found in thermal states [33], XY models [31] and in gapless systems 



The mutual entropy of this system, computed in (50) tells us that the work needed to 
erase all correlations [36] between two different blocks in the AKLT ground state decay 
exponentially to zero in the limit of infinite length. In this sense, all correlations between 
two blocks located infinitely far apart vanish (for non entangled boundary spins) which is 
expected from the thermodinamic limit of a gapped system. 

We also studied other different boundary conditions which, altogether with the results 
found in [37] , agree in the limit of infinite separation, except in the case when we start with 
an entangled pair at the spin 1 boundary of a free end AKLT ground state. In that case as 
shown in [37], the entanglement reduces to the entanglement of the Bell pair created from 
the spin 1/2 virtual particles which remain free after the antisymetrization between different 
neighbor sites. 

Our result is in agreement with the fact that ID AKLT chains alone are not sufficient 
for universal quantum computation. This is due to the vanishing negativity between two 
different non adjacent blocks. Still further coupling of many such chains can in principle 
implement quantum computation as shown in [2"3"] . 

Acknowledgements R. S. acknowledges the Fulbright-Conicyt Fellowship. V.K. acknowl- 
edges the Grant DMS-0905744. 



VI. APPENDIX 

A. Classical Variable representation 

A known representation of the boson algebra introduced in the first section is given by 
a\ = Ui, ai = b\ = Vi, hi = 7^7. In this representation the spin operators read 



d _ d I f d d \ 

St = U ^ S ^ Vi ^ S ^2{ Ui d^- V ^J (51) 

Due to the rotational invariance of the AKLT model, it's useful to choose [22] 
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Ul = e ^ 2 cos |, Vi = e-^' 2 sin |, (52) 

where 9 and <fi parametrize the unit sphere, with 9 e [0, ir] 9 = being the positive z axis 
and <f) e [0,2n]. 

The condition aja^ + fejfej = 2S imposes a restriction on the functions allowed to form 
spin states, namely 

Ui^- + v~) f(ui, Vi) = Sf(u h v t ) (53) 



2 V dm dvi 



The solution to (53) is f(u,v) = J2k fkU k v 2S ~ k , a polynomial of degree 25 in u and v, 
with a k an arbitrary constant. The inner product becomes (g\f) = J 1p-g(u,v)f(u,v) where 
Q is the solid angle over the sphere, g is the complex conjugate of g. In the subspace of 
degree 2S polynomials the matrix elements for S + are 



d f 

(g(u,v)\S + f(u,v)) = (g(u,v)\u—f(u,v)) = ^a,-a fc / 

V j,k J 



^_j_9S_» 9 h, no t 

— v?v ' b ] u—u V 

4:71 OV 



= ^2aja k 5 k+ld B(k + 2,2S - k){2S - k). (54) 

j,k 

with B(x,y) the beta function. Following [22], we introduce the classical variable represen- 
tation of S + as S£ = 2(S + l)uv. The matrix elements for this operator are 



(g(u,v)\S+f(u,v)) = (g(u 7 v)\2(S + l)uvf(u,v)) (55) 

dVL 

4tX 



j,k 



-u j v 2S - j+1 u h+1 v 2S - k 



= ^2a j a k S k+1J B(k + 2,2S -k + l)(2S + 2). (56) 

Now, writing the beta function in terms of the gamma function, and using T(z + 1) = 
zT(z), we have 
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then (g(u,v)\Stff(u,v)) = (g(u,v)\S + f(u,v)). The classical representation of the spin op- 
erators is S£ = (2S + 2)uv, S~ t = (2S + 2)vu, and — (S + \){uu — vv). This expressions 



provide the same matrix elements as the operators (51), as can be shown easily from the 
definitions. 

Now using the relations a = u = e ? ^ 2 cos|,a^ — u — e~* (?i//2 cos^ and b = v 



2' 

e -*0/2 g j n |^t — y — e l ^' 2 sin|, we can also prove that a similar relation holds for the 



overlap between states satisfying (53) 



(g(a,b)\f(a,b)) 



J§g(u,v)f(u,v) 



V(9\9)(f\f) y/j%\g(u,v)\*f%\f(u,v)\' 



(57) 



The state \i l+L l ) containing L sites fulfills the condition (53 ) at every lattice point, except 
at the boundary sites i and i + L — 1. The ground states |Af ) of the bulk Hamiltonian defined 



by \Afj) = TfJ(i,j + 1) introduced in satisfy the relation (53) at each lattice site. 

The norm of the VBS state Q is in this language 



(VBS|VBS) 



'N+l 

n 

,i=0 



dVLi 

47T 



N 



i=0 



1. 



(58) 



with Cli being the radial vector on the unit sphere Cli = (sin 9^ cos 0,, sin 9i sin 0,, cos 9j). 
The norm of the states \A^) = T^(i,i + L — 1) |/ +L_1 ) composed by L sites is then 



(AJA,) = jfl ^(i - Ck ■ <W)V(i, L)%(i, l) = \ ^ 



l + s L 



Sau, (59) 



with s. 



-1,-1, 3,-1). Here we have introduced 7^, the classical analog to the operators 



T„, defined as 



= <Pi(<Tn)ab<p b j, With tp\ = U i: (fit = Vi. 



(60) 



B. Identities 



All the identities that we have use in this work can be obtained from the basic identity, 
(repeated indices are summed) 
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= -\{-irV^k){a,) ab , (6i) 

which can be checked directly by inspection, and makes use of the fact that the matrices 
cr M form a basis of the GL(2, C) group. Then for two boundary operators (these operators 
appear naturally in the boundary between two different blocks in the bulk) d\ = T 2 \i, i + 
and = T 2 \j,j + 1), we have 



d\d) = ^(a 2 WliVf(^)c^+i 



now applying the identity (61) twice and using the fact that (<7 u ) a d = [&v)da which is also 



equal to -{-\) v g vv ' (o v ,) da , we get 



djdj = I(-l)/y-V(* + !> 3)Tj(i, 3 + l)Tr(<r 2 <v7 2 <v)- 

The last term in the above expression is the trace of 4 matrices. To compute it we can use 
that 0"2<7 M cr2 = (—1)^0"^ and that T^cx^ov) = 2g^ u obtaining 

djd} = -^(1 + 1,^(1, j + 1) (62) 

Using this identities it is possible to generate all the identities for any number of boundary 
operators d. In this work, we used the identity for three and four d operators 



>i»r k ~= -\{-l) v g vv 'V{i + l,j)Tj(j + l.k)Tj(i, k + 1) x Tr(a^a 2 a x ) 

(63) 

dldjdldj = t^lg^'T^i + l,j)T u \j + 1, k)T p \k + l,l)T x \i,l + 1) x Tr(a,a„,a p a x ). 

(64) 

here o p = (i, a%, cr 2 , 03) where 0*. are the three Pauli matrices. We also define a p = 
(— z, 01, (72, 03)- To compute the traces of the Pauli Matrices, we use the following tricks. 



<W + cr u a p = 25 n„. 
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The a object is a generator of the Lorentz transformations in Euclidean space, so it 
satisfies the Euclidean Lorentz algebra 

[(Tpv, a <xp] = 1{&va<J ii$ - <^/3<V* + Spp&va ~ ^q^^ (65) 

Further identities can be derived using the Dirac matrices technology, namely, in Eu- 
clidean space we have 









<i o\ 




|. 7 5 =| 




- 

U(Th 


° 1 




10 -I 



7m = h 75 = (66) 

J \0 -I J 

then we can compute the trace of four Pauli matrices using the identities for Dirac matrices 
[28] . and projecting out the lower block with the chiral projector (1 — 7s)/2. We have for 
example 



Tr(cr^CT V <7 P <7\) = ^Tr(7 M 7 z ,7 P 7 A (l - 7 5 )) = 2(5^ P \ + <V<V - Sn P S v \ + e^px) 



C. General Results 

1. Block separation L > 1 



In section |III| we found an explicit expression for the density matrix of two blocks of 
length L\ and L 2 separated by L sites. In that section we presented the asymptotic results 
for the eigenvalues of pab- The result for any L\, L 2 > 1 is given in terms of the following 
quantities 



11. _ / 1\ L \ . if, ( l x 1 



X ^ = l[ 1 + 3 {-s) J> A 1 (L) = i (l 

From those quantities we define A o = Ao(£i)A (L 2 ), An = Ai(Li)Ai(L 2 ) and A10 = 
Ao(L 1 )A 1 (L 2 ) + A (L 2 )A 1 (L 1 ). 

The characteristic polynomial associated to the density matrix pab, p(Y) = det(pAB — Y) 
is p(Y) = pi(Y) 5 p2(Y)p3(Y) 3 , where Pk(Y) is a polynomial of degree k on K, given by 

(z = (-3)- L ) 
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p 1 (Y) = Y-(l-z)Xu, 

p 2 {Y) = Y 2 - (A 00 + (1 + 2z)X n )Y + (1 - z)(l + 3z)A 00 An, 

p 3 (Y) = Y 3 - (A 10 + A n (l + z))Y 2 + [(1 + z)A 00 + (1 + 2z)X w ]{\ - z)X xl Y 

- (1 - z) 2 {\ + 3^)AooAn = Y 3 + bY 2 + cY + d. (67) 

We define q = ^(2& 3 — 96c + 27c?), p = |(3c — b 2 ). The eigenvalues of pab are the 
solutions to P(y) = 0. They are 



y = (1 — z)Xn five-fold degeneracy, (68) 
V = \ ((Aoo + (1 + 2z)X xl ) ± ^(Aoo + (1 + 2z)A n ) 2 - 4(1 - z)(l + 3z)A 00 Aii) , (69) 

y = 2i/-f cos ( J arccos f -- j + ^ j - J, (A; = 0, 1, 2) triple deg. (70) 



3 \3 \2p\ pj 3)3 

2. Adjacent blocks 



In section III, we computed the transposed density matrix of a system consisting of two 
blocks inside the VBS state, A and B, of length L\ and L 2 respectively. The spins which do 
not belong to AUB have been traced away. In the general case when the blocks are separated 
by L sites we could prove that the negativity vanishes for L > 1, being the only nontrivial 
case when L = 0. In that case, the negativity in the asymptotic limit L\ — > oo, L 2 — > oo is 
given by 

^-=5-i((-ir + (-ir)- (n) 

The decay in the thermodynamic limit is twice as fast compared to the usual decay of 
the spin correlations, a feature that is already seen in the case of the negativity of the pure 
system studied before. 

The logarithmic negativity also show this behavior, for L 1; L 2 ^> 1 

1 A / 1 x 2 ^ 2 



^ =1 -iE(2jH~) ' (72) 
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In the of adjacent blocks, for any L%, L 2 > 1, the characteristic polynomial associated to 
the transposed density matrix p^, p(y) = det(p^g — Iy) is p(y) = Pi{y) 5 p2{y)p3{y) 3 > where 
pk(y) is a polynomial of degree k on y, given by (defined in terms of A o, Ai , An) 

p x {y) =j/-2Au(l,2), 

ft(y) = ?/ 2 - y(Aoo(l, 2) - A n (l, 2)) - 4A 00 (1, 2)A n (l, 2), 

ps(y) = y 3 - A 10 (l, 2)y 2 - 2A 10 (1, 2)A n (l, 2)y + 8A 00 (1, 2)A n (l, 2) 2 = y 3 + % 2 + cy + d. 

(73) 

this polynomials are related with the polynomials of previous section by taking z = — 1 in 



eq. (67). Out of the 16 eigenvalues (y n ), 4 are negative, with y\ (no degeneracy) given by 



the expression 



l - (a 00 (1, 2) - A n (l, 2) - ^00(1, 2) 2 + 14A 00 (1, 2)A n (l, 2) + A n (l, 2)^ (74) 
and y2 (triple degeneracy) given by 



Vi 2 



y2 = " 2 V "3 sin G arccos (I V -p) + 1) " 5 (75) 

The negativity of the system is then A/" = — (yi + 3?/2), while the logarithmic negativity 
is given by Ej^ = log 2 (l — 2(yi + 3^))- In the special case when L\ = L 2 = I, the negativity 
simplifies to (using x = (— 3)~ z ) 

Af(l) = + x 2 - + Ax + 2x 2 - 4x 3 + 13a; 4 - ^ ^{l + 3x)(l - x)^j , (76) 

1 3 

~ - - -(x 2 - x 3 ), for sc < 1. (77) 
The logarithmic negativity is given by Ej^(pab) — 1 — 2 in(2) ~ 
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